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Abstract: While amongst the most luminous objects in the universe, many details regarding the 
inner structure of quasars remain unknown. One such area is the mechanism promoting increased 
polarisation in the broad absorption line troughs of certain quasars. This study shows how microlensing 
can be used to differentiate between two popular models that explain such polarisation through a 
realistic computational analysis. By examining a statistical ensemble of correlation data between two 
observables (namely image brightness and polarisation of the flux coming from the quasar), it was 
found that through spectropolarimetric monitoring it would be possible to discern between a model 
with an external scattering region and a model without one. 
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1 Introduction 



gravitational lensing — polarisation — quasars: individual 



Sections 4 and 5 respectively. 
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Quasars are amongst the most luminous objects in the 
universe, radiating most of their energy from within 
a small continuum emitting region only ~ 1 pc in ex- 
tent. Such a source at cosmological distances subtends 
the order of microarcseconds, well below the highest 
angular resolution obtainable with modern telescopes. 
However, gravitational microlensing can be employed 
to reveal the structure at the heart of quasars, with 
magnification due to individual stars revealing infor- 
mation ab out the size of the continuum emitting accre- 
tion disk {W ambsganss, Sch neider, Paczynskilll990l: 



2.1 



Background 
BAL Quasars and Polarisation 



Approximately 10-20% of optically selected quasars 
have been found to exhibit broad absorption troughs in 
resona nt lines blue- ward of t he corresponding em ission 
lines (jHewett fe Foltzll2003l : iReichard et alJl2003h . ex- 
hibiting bulk outflow ve locities of around 5000-30000 
km s~ 1 (lTurnshekl [i984 N ). These BAL quasars 



Witt fc Ma o 1994: Hawki ns fe Tavlorlll997l:lGould & Mirald&lfemlde]o consist of only radio- quiet or radio- intermediate 



19971) and broad emission line region (IWvithe &; Loebl 
20021 : lAbaias et aT1l2002l : iLewis fe Ibatall2004T ). 

While the microlensing advances have improved 
our understanding of quasar structure, the picture is 
far from complete, especially with regards to areas 
such as the mechanism promoting jet activity [ which in 
turn d i ctates rad i o- loud n ess or -softne s s, e.g . iKuncid 
(jl999f ): ICattaneol (|2002T ): iFender et al.1 d2Q04h1 and its 
relation to prominent absor bing regions fe.g. lMurrav et a"D 
119951 : ILewis fe Bellel Il998h . In particular, the exact 
mechanism which promotes increased polarisation in 
the broad absorption line (BAL) troughs of quasar 
spectra is not well understood . This paper will ex- 
tend recent studies ([Lewis &; Be lle 1998; Belle & Lewis 
2000) to computationally simulate how gravitational 
microlensing can be used to discern between the two 
main models for polarisation enhancement and hence 
probe the scales of structure in the absorbing/scattering 
regions. Section 2 discusses of the details of BAL 
quasars, microlensing, and the quadruply imaged quasar 
H1413+1143, the system which is the ideal observa- 
tional candidate for this study. Section 3 details the 
approach to the simulations, including the construc- 
tion of the source profiles and polarisation maps, with 
the results and conclusions of this study presented in 



sources, but the discovery of a radio-loud BAL quasar 
suggests the phenomenon occur s throughout the quasar 
population ([Becker et al.lll997h . Current theories sug- 
gest an orientation-based unification schem e to explain 
the oc currence of the BAL in quasars (e.g. lAntonuccil 
(1993)), and of interest to this study is the scatter- 
ing structure which increases polarisation within BAL 
troughs. 

Figure [T] shows the two competing models for ex- 
plaining the enhanced polarisation in the broad ab- 

sorpt ion lines ([Cohen et al . 1995; Ogle 1997; Schmidt & Hines 
1999). For Model A, radiation from the nuclear region 
only travels along Path A through the BAL clouds. 
This radiation is then enhanced by resonant scatter- 
ing into the line of sight, increasing the polarisation 
within the absorption troughs. For Model B, radia- 
tion travels along Paths A and B. However, in this 
case Path A does not necessarily introduce any in- 
crease in polarisation. Radiation travelling along Path 
B is scattered into the line o f sight, most likely from a 

region of electrons and dust (lAntonucci &; Millerll 19851: 

[ Goodrich fc Millerl[l995l : lGallagher et al .111 9991 : iBrandt et all 
I1999T ). Here, the increased polarisation in the absorp- 
tion troughs is due to the reduced amount of unpo- 
larised flux coming directly from the continuum along 
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Path A (the polarised flux remains relatively free of 
absorption as it avoids the BAL clouds). 

2.2 Microlensing 

The details of gravitatio nal lensing can be fo und in re- 
cent review articles (e.g. Wambsganss 1998) and only 
the salient features will be discussed here. Gravita- 
tional lensing occurs when light rays from a distant 
source pass near a massive object and suffer achro- 
matic deflection. The deflection angle a of such a light 
ray passing at a distance r from an object of mass m 
is given by Equation [T]- i.e. light rays will follow null 
geodesies in the presence of massive objects. 
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The measurement of this deflection was one of the first 
key observational tests of Einstein's general relativity. 

When considering the deflections due to individual 
galaxies, it is seen that multiple light paths can con- 
nect a source with an observer, resulting in multiple 
imaging on the scale of arcseconds. However, when 
small-scale granularity in the distribution of galactic 
matter (stars, planets, black holes) is considered, it is 
seen that these macrolensed images are actually com- 
posed of a myriad of unresolvable (~ 10 -6 arcsec) 
microimages due to im aging by stellar- mass objects 
(jChang Refsdal 1979). While these microimages are 
unresolvable, the stellar mass objects can introduce 
magnification of a background source, and the motions 
of the lensing stars can produce significant fluctuations 
into the observed brightness of the macroimage. The 
characteristic scale length of microlensing is the Ein- 
stein radius (ER). For a point mass lensing, and a per- 
fect alignment of source (s), lens (1) and observer (o), 
the result would be a circular image known as an Ein- 
stein rin^Q at the Einstein Radius (ER). The physical 
projection of this radius onto the source is given by 
Equation El where D is the angular diameter distance. 



ER 
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the mean of these two values, z = 1.55, to represent 
the re dshift of the lensing g alaxy. One of the images, 
D, fsee llChartas et all (|2004f )l. has been seen to exhibit 
variability consistent with microlensing (I Angonin et al.l 
ll99d:|Kavser et alJl990t l0stensenl 19971 ). with lHutsemekersI 
(|1993h suggesting that prominent differences in the ab- 
sorption profiles of the images might be due to selective 
microlensing of absorbing clouds, with the scale size of 
these clouds being smaller than the continuum- forming 
region. It was conceded, however, that this would re- 
quire a very precise lensing configuration. Further- 
more, the polarisation for the summed images has been 

seen to vary in the blue wing of the Cj y A 1549 emission 

line, with fluctuati ons between ~10% (lLamv &; Hutsemekersl 
12004 ) and -20% flSchmidt fe Hinesl Il999h . Interest- 
ingly, these studies also found that continuum polari- 
sation near the Cjy feature was —2%, indicating some 
polarisat ion is still occurring away from the absorption 
troughs (jWang. Wang, fe Wandl2005f ). 



There seems to be growing support in the litera- 
ture that an external, asymmetric scattering region is 
responsible for this polarisation increase + (Model B 
in Section [27TT). A recent HST study of H1413+1143 
(|Chae et al.ll200lh has indicated that the size scale for 
such a region lies approximately between the Einstein 
ring size (i.e. Einstein diameter in the source plane) 
and 10 18 L46 5 cm (where L46 is the lensed quasar lu- 
minosity in units of 10 46 ergss -1 ). The minimum 
value comes from the requirement that the scatter- 
ing region should not be as susceptible to microlensing 
as the central nucleus, whereas the maximum value 
is simply an estimate for th e size scale of the broad 
emission line region (BELR) (Mur ray fc Cha ing 1998; 
iKaspi et al.ll2000h . If the scattering region were any 
larger then the differences seen between macrolensed 
images would not be as great. Furtherm ore, recent 
data from Chandra ([Chartas et al.l I2004T ) appears to 
support this model. 



By applying more realistic models for both the 
BAL region of H1413+1143 and the distribution of 



The importance of this length scale is that objects 
smaller than it are much more susceptible to large mag- 
nificat ions ,wlril£_ojDje^ — ■ - ..^^ ( .... 

nification (IWambsganss k, Paczvnskilll99ll : ISchneider. Ehler &t&sBgil^ie lensin g galaxy, the observe d spectral varia- 
1999; IWambs ganss 1992^ tions examined bv lHutsemekersI (|l993h could be readily 

reproduced, removing the need f or precise lensing con- 
figurations (Eewis & Belle 1998). This idea was fur- 
thered by Belle & Lewis ([20001 ) who investigated the 
role that Models A and B (see Figure [TJ might play 
in explaining the increased polarisation within BAL 
troughs. However, this previous study only tested one 
configuration with a scattering region scale size of half 
an Einstein radius. Hence, how secure are the conclu- 



2.3 H1413+1143 

H1413+1143 (the Cloverleaf) consists of 4 images of 
a z — 2.55 quasar with angular separations of 0.77 t o 
1.36 arcsec (|Magain et al.lll988l : lTurnshek et al.lll997f ). 
The lensing galaxy has only recently been identified 
dKneib. Alloin. fc Pello 1998; C hantry fe MagairfeOO/t) . 
with iMagain et al.l (|l988h and lAngonin et all (jl99Ch 
identifying two prominent absorption systems at z — 
1.438 and 1.661; for the purposes o f this study, and 
for consistency with Lewis & Belle (1998), we adopt 



1 Despite being first proposed bv iChwolsonl (|192 



sions drawn by IChae et al.l (2001) and Chartas et al. 
(|2004r ). and could Model A produce polarisation vari- 
ations that could be misinterpreted as Model B? To 
this end, the remainder of this paper examines detailed 
simulations of both models and the predictions they 
make for polarisation variability. 
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Figure 1: The two potential m odels for the scattering and absorption structure (see explanation in Section 
O also lBelle fc Lewis! (|2QQQh ), 



3 Method 
3.1 Ray Tracing 

When considering microlensing at cosmological scales, 
many stars influence the path of light through a galaxy, 
and the single, isolated lens approximation which works 
well in Galactic microlensing must be abandoned. Such 
a situation is analytically intractable and numerical 
techniques must be employed. This study employs the 
inverse /backwards ray- shooting technique (jKavser et al.l 
1990; Wambsganss 1990) which involves 'shooting' light 
rays from the observer to the lens plane, calculating 
the deflection due to individual stars in the lensing 
galaxy, and then collecting these rays in the pixels of 
the source plane where they eventually hit, forming a 
magnification map. 

There are two main parameters required to model 
the microlensing which are dependent upon the mass 
distribution in the lensing galaxy - the dimensionless 
surface mass density a (or optical depth) and the shear 
7 due to the large scale matter distribution. As these 
parameters are not strongly constrained for the Clover- 
leaf, this study has investigated 4 different models: 
cr = 7 = 0.25, 0.4, 0.6 and 0.75. These were chosen be- 
cause they represent a singular isothermal sphere, rea- 
sonably approximating the rang e of possible mass dis- 
tribu tions for a lensing galaxy (Schneider, Ehle rs. &; Falcol 
1999). Figure [2] shows example magnification maps for 
regions derived using the ray tracing method, with side 
lengths of 10 ER for the 4 cases of a and 7 above. 

These maps can be scaled to physical distance us- 
ing the Einstein radius (Equation [2}. Assuming a stan- 



dard cos mological model with h = 0.73, A = 0.76 and 
Q = 0.24 ([Spergel et al.l l2006). and the lens and source 
redshifts given in Section ED ER^2.73xl0 16 cm. 

The magnifications due to microlensing fluctuate 
about a mean theoretical value given by 

/ x th =[(l- ff a )- 7 a ]- 1 (3) 

which is the magnification an image would suffer if 
the macrolens was composed of solely smooth matter. 
Hence, generally, the mean magnification over a large 
enough area should tend to this theoretically expected 
value. However, due to statistical variance, the mean 
value within an individual magnification map which is 
relatively small can deviate from this expected value 
(i.e. the smaller the region chosen, the larger the pos- 
sible deviation from the mean theoretical magnifica- 
tion) . 

3.2 Defining the Source 

In order to determine the effect of microlensing on 
each model from Section [2.31 it was necessary to con- 
struct the image that the magnification map would 
'see'. This is shown schematically in Figure [3] and dis- 
cussed in detail below. 

3.2.1 Absorption and Polarisation 

In order to examine a BAL trough, Equation [4] was 
used to simulate the blue wing of Cjy A 1549, where 
A cc = 1505 A is the wavelength at which 50% of the 
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Figure 2: Example magnification maps for <j = 7 = 0.25 (top left), 0.4 (top right), 0.6 (bottom left), and 
0.75 (bottom right). Darker regions indicate higher magnification, with the sharp boundaries denoting 
caustics of the map. Each have dimension 10x10 ER. 




Figure 3: Schematic view of the BAL quasar continuum and scattering region as seen by the microlensing 
magnification for both scattering scenarios. This is based on Figure Q] and is discussed further in Sections 
[3221 and [3231 
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light from the continuum core is absorbed (correspond- 
ing to a bulk outflow at 50% absorption of ^8500 
krns" 1 ) and 2AA = 8.75 A is the 'width' of the ab- 
sorption feature betwe en 27-73% absorption (|Turn shek 
fl995l:IWevmanr][T995l). 



A (A) 



1 + exp [(A - A cc ) /AX] 



(4) 



With regard to polarisation, as discussed in the fol- 
lowing sections, the polarisation at any pixel was di- 
rectly related to the absorption at that pixel. This 
could take on any value in the range from 0-20%. This 
would tend to generate conservative results for this 
study (note that, as mentioned in Section [2731 the po- 
larisation of the entire composite image could be as 
high as 20%); however, any general trends established 
here would still apply if a larger range for polarisation 
was used. 



3.2.2 Model A: Scattering Within the BAL 
Region 

For this study the quasar continuum source was mod- 
elled as a two-dimensional Gaussian surface with a 
brigh tness radius of 10 15 cm (Rees 1984; B lanford. Netzer. 
1990) on a 64x64 pixel grid. The extent of this grid 
was defined so that the value of the Gaussian at the 
edges would be 10% or less than the value at the peak, 
giving a side length of 3.65xl0 15 cm (this kept the size 
of the grid down while allowing for good resolution). 
Hence, this Gaussian represents the unabsorbed con- 
tinuum flux and in order to represent absorption, an 
absorption matrix (discussed shortly) would then be 
(dot) multiplied with the source grid (i.e. this absorp- 
tion matrix is the same size as the source matrix). 

In order to represent the scale length of the inho- 
mogeneities in the BAL region (this will be referred to 
as cloud size), cloud sizes were set as 3.65xl0 15 x4 _n 
cm, n=l,2,3. In this way, the smallest clouds were 
5.70xl0 13 cm in extent (n=3) while the largest clouds 
were 9.11xl0 14 cm (n=l). This was implemented by 
limiting the initial matrix size for the absorption ma- 
trix generator and then drawing this matrix out to fill 
a 64x64 grid . Note, this method was similar to that 
employed bv lLewis &; Belld (|l998h . 

The absorption matrix represents the degree of vari- 
ation of absorption about a particular wavelength (i.e. 
due to inhomogeneities in the density of the absorption 
region) . The distribution of values were represented by 
a Gaussian of specified width co centred at A C (X). A 
modification algorithm was developed to truncate this 
distribution outside 0-100%, and to reassign new val- 
ues for the truncated elements such that the overall 
distribution of matrix values had mean (A (A)) and 
width uj. The mean of the absorption distribution 
(A (A)) was chosen to be a function of wavelength as 
given by Equation [4] The algorithm was also designed 
so that if, for a particular (A (A)) and a;, A C (X) was lo- 
cated outside the 0-100% boundary, then A C (X) would 
be fixed at the boundary and u) varied so as to ensure 
that (A (A)) was the value specified. This allows, for 
example, for a distribution with (mean) absorption of 



10% and width of 50% to be generated so that all ma- 
trix values lie within 0-100%. For this paper, widths of 
5% and 50% were investigated for absorption at 25%, 
50% and 75%. Figure |4] shows how different cloud sizes 
were modelled with a particular set of absorption ma- 
trix parameters. 

In order to implement polarisation, this image was 
then split into two separate images according to Equa- 
tion [5] In this way, if Equation [6] was calculated at 
each pixel then the polarisation would be dependent 
on the particular absorption matrix value at that pixel. 
Here, Pix n are the values at corresponding pixels in 
each new image, Pol una b = 0.02, Polu m = 0.18, Abs is 
the absorption matrix value at that pixel, and Pol is 
the polarisation difference between the two images. In 
this way, if Abs=0% then Pol=2%, if A6s=50% then 
Pol=ll%, and if Abs=100% then Po/=20%. 



Pixi = O.bPtot (1 + PoLnab + Poli im x Abs) 

PlX 2 = O.bPtot (1 - Polunab ~ Poliim X Abs) 

Ptot = Pixi + Pix 2 

Pdif = Pix\ - Pix 2 (5) 



fc Woltierl 



Pol ■ 



Pdif 

Ptot 



(6) 



3.2.3 Model B: External Scattering Region 

This model was somewhat constrained by the maxi- 
mum magnification map that could be handled in terms 
of computer memory, so by setting this map with a side 
length of 100 ER (1024x1024 pixels) the pixel scale size 
was automatically set for any source models created. 
Fortunately this did not seriously affect the range of 
scattering region scale sizes that needed to be investi- 
gated. Unlike Model A, an absorption matrix was not 
required here (the continuum region is on the order of 
a few pixels), but rather absorption simply involved 
reducing the flux of the continuum while keeping the 
flux of the scattering region constant. 

This left two main variables to manipulate besides 
absorption: the scale size of the scattering region and 
the separation distance of this region from the contin- 
uum. Note that the latter is in fact intrinsically linked 
to the scale size of the scattering region in that at scale 
sizes approaching the upper limit, the separation dis- 
tance must be of the order of the scattering region ra- 
dius (geometries where the separation distance is less 
than the scattering region radius would give similar 
correlation results). The scattering region radii and 
separation distances investigated for this model are 
shown in Figure [5] This involved setting the source 
as a 64x64 grid for the smallest scattering region and 
200x200 for the larger two regions (resulting in much 
higher computation times). Also note that, like the 
continuum region in Model A, both the continuum and 
scattering regions in Model B were modelled using two- 
dimensional Gaussian surfaces. 

The method for achieving polarisation in this model 
is quite different to that of Model A. Here, polarisation 
is due to the reduced flux coming from the continuum 
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Figure 4: Model A. Left: Source with cloudsize 9.11xl0 14 cm (n=l), absorption at 25% and uo at 50%. 
Middle and Right: Same parameters but with cloudsize 2.28xl0 14 cm (n=2) and 5.70xl0 13 cm (n=3). 
Each have side length 3.65xl0 15 cm. 




Figure 5: This shows the 5 different sources considered as part of Model B at high absorption percentages 
(so that the continuum doesn't drown out the scattering region in the images). Scattering region radii are 
3.00xl0 16 cm with separation distances of 6.03xl0 16 (top left) and 1.66xl0 17 cm (top right), 7.70xl0 16 cm 
with separation distances of 2.00xl0 17 (bottom left) and 5.39xl0 17 cm (bottom middle), and 2.00xl0 17 
cm with a separation distance of 3.70xl0 17 cm (bottom right). Images in the top row have side length 
1.71xl0 17 cm whereas images in the bottom row have side length 5.34xl0 17 cm. Note the continuum in 
the bottom left corners. 
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so that, unlike Model A, a linear increase in absorption 
will not create a linear increase in polarisation. For 
this reason, absorption values of 77%, 91% and 97% 
were chosen so as to represent the same mean polari- 
sation values used in Model A. To generate this source, 
the continuum region flux is diminished in accordance 
with the absorption value and then split between two 
images. The scattering region is then added to each 
image so that they satisfy Equation [6] for a constant 
polarisation of 20%. In addition, the relative flux of 
the two sources can then be set by the requirement 
that the polarisation at 0% absorption is 2% (this then 
increases non-linearly to 20% polarisation at 100% ab- 
sorption) . 



3.3 Convolution 



In order to determine the magnification distributions 
of the various source images, these images need to be 
convolved with the magnification maps. In doing this 
it was important to ensure that the physical scale cor- 
responding to the pixel length of both source and mag- 
nification map was identical. For Model A this was de- 
termined by the design of the source, requiring a mag- 
nification map of 2.136x2.136 ER (1024x1024 pixels). 
As this region is reasonably small it was important to 
recognise that a number of magnification maps would 
need to be used to obtain data (see Section l3J"|) . Ne- 
glecting this could reveal incorrect trends in any cor- 
relation data (in particular it would reduce the range 
of magnifications significantly). For this study, two 
maps of varying mean magnification were used to en- 
sure that possible trends at both large and small mag- 
nification would not be overlooked (e.g. a trend at 
low magnification might be overlooked if dealing with 
a magnification map excessively covered with intense 
caustic structures); this is discussed further in Section 
14.11 Note that this problem does not exist for Model 
B because the map is large enough that, on the whole, 
it is reasonably homogeneous and so has a mean that 
approaches the theoretical expectation value. As men- 
tioned previously, the process for Model B was the 
reverse of that for Model A, whereby the largest mag- 
nification map of 100x100 ER was created in order for 
the source to then integrate with it. 

In order to determine the polarisation following 
convolution (using the two polarised images formed 
for Model A or B) Equation [6] was used. These values 
could then be compared directly with the magnifica- 
tion values obtained from convolving the magnification 
map with an unabsorbed, unpolarised version of the 
corresponding source (these values were divided by the 
mean flux of the original source in order to obtain the 
correct magnification values). Correlation maps could 
then be obtained. Finally, note that for this paper all 
convolutions were computed without any zero-padded 
edges, so that areas contaminated by the overlaps were 
neglected. 



4 Results 

Before looking at any correlation data it is insightful to 
look at the light curves to see just how the mean mag- 
nification and polarisation vary, as shown in Figures 
[6] and [71 These light curves are generated by taking 
a slice through the convolved data, representing what 
an observer would see as caustics pass over the source. 

The time scales shown in these figures have been 
calculated using Equation [71 where the scale size of 
the source is /15IO 15 cm, zi is the redshift of the lens- 
ing galaxy (z ~ 1.55), and 300^300 kms _1 is the ve- 
locity of the microlen sing stars across the line of sight 
(jLewis fe Belldll998h . 



Di 1 



1 + zi D s V300 



(7) 



However, it is worth noting that the uncertainty 
in the redshift of the lensing galaxy results in an un- 
certainty on the Einstein radius in the source plane 
and hence the physical size of the inferred absorption 
and scattering regions. In considering potential lens 
redshifts between z — 1 and z — 2 (assuming the cos- 
mology from Section l3.1[) it can be shown that the size 
of the Einstein radius varies roughly linearly between 
+30% and -30% of the value at z ~ 1.55. In other 
words the time scales shown in Figures [6] and [7] would 
vary by these same percentages. 

These results hint at some interesting correlations 
between magnification and polarisation. In Model A, 
the polarisation fluctuations tend to occur at caustic 
crossings and the polarisation seems to oscillate about 
a mean value. In Model B it is seen that decreases 
in polarisation tend to occur at higher magnification 
and without any significant oscillation. This can be 
seen in the top image of Figure [71 whereby the polar- 
isation during each microlensing event decreases sig- 
nificantly rather than oscillating about any particular 
value. With regard to possible observational testing 
this last point is of crucial significance because it indi- 
cates that as few as one major microlensing event may 
be used to differentiate between the two models. It is 
also worth noting that the same number of microlens- 
ing magnification events can be statistically expected 
for both models within a 20 year period, with this be- 
ing related to the time taken for caustic structures to 
sweep the continuum source. As we will show in the 
next section, these trends are confirmed in an analysis 
of the correlations between magnification and polari- 
sation for a large sample of microlensing scenarios. 

4.1 Trends in Model A 

The correlations between polarisation and magnifica- 
tion for the scenarios described in Section 13.2.21 are 
shown in Figure [8] each panel presents a grey scale 
map of relative probability, with the continuum mag- 
nification on the x-axis and the degree of polarisation 
on the y-axis. Hence this shows trends between any 
fluctuations in the two. Before discussing the trends 
seen it should be noted that the apparent decrease in 
mean polarisation with increased absorption width (in 
this model) is in fact an artifact of the computational 
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Figure 6: Model A. Left: Convolution of unabsorbed source with 2.136x2.136 ER magnification map 
(a = 7 = 0.4). Middle: Resultant polarisation after application of Equation [6] to two convolutions of 
the polarised source (absorption=50%, width=50%). Both have side length 5.47xl0 16 cm. Right: Light 
curves for magnification and polarisation, as indicated by line from top to bottom in previous images. 




Figure 7: Model B. Top Left: Convolution of unabsorbed source with 100x100 ER magnification map 
(scattering region radius = 3xl0 16 cm, separation distance = 1.66xl0 17 cm and cr = 7 = 0.4). Top 
Middle: Resultant polarisation after application of Equation [6] to two convolutions of the polarised source 
(absorption=91%). Both have side length 2.56xl0 18 cm. Top Right: Zoomed in view of bottom plot, 
corresponding to first 100 years. Bottom: Light curves for magnification and polarisation, as indicated 
by line from top to bottom in top left and top middle images. Note that the axes are slightly different to 
those in Figure [U 
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model. This artifact does not influence the overall 
trends obtained. 

Firstly, it is apparent that the choice of matter dis- 
tribution parameters for the lensing galaxy (i.e. a and 
7) do not significantly affect the trends seen in these 
probability distributions (especially with the knowl- 
edge that any real data would be both incomplete and 
have error associated with it - only the most general 
of trends would be noticeable) . The same can be gen- 
erally said for the choice of cloud size, although it can 
be seen that the range of polarisations observed at 
any particular value of magnification decreases with 
decreasing cloud size. It is also worth noting that 
the choice of mean absorption does not induce any 
trends into the data apart from an obvious increase in 
mean polarisation, indicating that possible spectropo- 
larimetric monitoring would not reveal different trends 
at different wavelengths within the Cjy trough (recall 
Equation |4}. In addition, by using two sets of data for 
each correlation map (corresponding to a small and 
large mean for the magnification map used, see Sec- 
tion [33}, as seen more noticeably for the cases where 
n=l and cj=50%, it can be seen that the overall trends 
developed do not alter with lower or higher magnifica- 
tion. 

All of the scenarios tested seem to support the gen- 
eral trend that an increase in magnification of the con- 
tinuum flux (i.e. the unabsorbed continuum) will not 
indicate an increase or decrease in mean polarisation, 
but rather that the polarisation will remain reasonably 
constant throughout any microlensing activity. For all 
cases it is clear that increasing the width of absorption 
will increase the range of polarisations found at any 
particular magnification (the largest ranges are found 
around a magnification of 1, which of course is the 
approximate mean value of the magnification map). 
This then implies that if the source was observed over 
several months, the largest range of polarisations mea- 
sured would have been occurring when the mean mag- 
nification of the continuum flux was around the theo- 
retical mean. In other words, provided this model rep- 
resents the view of the quasar, then the largest fluctua- 
tions in polarisation would occur when the continuum 
flux was at its time-averaged mean. 

4.2 Trends in Model B 

The correlations between polarisation and magnifica- 
tion for the scenarios described in Section 13.2.31 are 
shown in Figure [9] It is immediately apparent that 
the trends seen in this model are very different from 
those seen in Model A, and at low magnification (for 
all scenarios tested), the polarisation values are much 
greater than at higher magnifications. Not surpris- 
ingly this is because of the the geometry of the source 
as only the continuum source is strongly magnified, 
whereas the scattering region is effectively immune to 
the influence of the microlenses. 

As with Model A, differences between mass distri- 
butions parameters for the lensing galaxy do not signif- 
icantly affect the overall trends seen in the data. How- 
ever, unlike Model A it can be seen that the choice of 
mean absorption does introduce a specific trend (other 



than simply increasing the mean polarisation) ; namely 
that at low absorption the mean polarisation drops 
with increasing magnification while at higher absorp- 
tion the mean polarisation tends not to drop as much. 
In fact, at very high absorption this correlation be- 
gins to look like a horizontal line. This is illustrated 
further in Figure [10] (as per a row of Figure [9} for 
a = 7 = 0.4 (this trend is the same for the other mass 
distributions), where absorption is now 99.7% (corre- 
sponding to a wavelength of ^1530 A). This result is 
important because it indicates that spectropolarimet- 
ric monitoring would reveal different trends at different 
wavelengths within the Cjy trough. 

Finally, it can be seen that by increasing the size of 
the scattering region the relationship between magnifi- 
cation and polarisation becomes more one-to-one (less 
broad). This occurs because, as the ratio between con- 
tinuum and scattering region scale sizes decreases, the 
unpolarised continuum is allowed to be magnified in a 
more dominant manner thus reducing the overall po- 
larisation. It can also be seen that variation in the 
separation distance between continuum and scattering 
region centre does not introduce any significant trends 
into the data. This makes sense, because as soon as 
these two regions are asymmetrically separated by at 
least the width of a caustic then the manner in which 
they are magnified becomes independent, yielding the 
similar results seen here. 



5 Conclusions 

This study has investigated the role that microlensing 
can play in differentiating between two popular models 
of how quasar BAL troughs are polarised. Using the 
macrolensed and microlensed quasar H1413+1143 as a 
case study, two computational models of polarisation 
were developed in order to investigate how magnifica- 
tion and polarisation variation detected by an observer 
could be used to differentiate between them. The re- 
sults showed that the correlations between these two 
observables would be easily discernible between mod- 
els. Two main differences were identified; for Model A, 
during a single microlensing magnification event, the 
polarisation at any wavelength within the Cjy trough 
was seen to oscillate about a mean value, while away 
from high magnification events the polarisation was 
found to be relatively constant. This was found to 
be in stark contrast with Model B, whereby the mean 
polarisation at any particular wavelength within the 
trough was found to rise and fall in anti-correlated 
fashion with magnification. The second main differ- 
ence was that if polarisation and magnification varia- 
tion were monitored at various wavelengths within the 
Cjy trough, then significantly different trends would 
be seen in this data between wavelengths for Model 
B but not Model A. Hence spectropolarimetric moni- 
toring through even a single high magnification event 
would provide constraints on the underlying scatter- 
ing geometry. Given that caustic crossing times for 
H1413+1143 would be of order months to years, an 
observational campaign monitoring this system on a 
weekly to monthly basis would be required to differen- 
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Figure 8: Correlations for Model A. Note that the scale reflects relative probability as denoted by the 
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(this has been normalised with respect to the theoretical mean magnification n t h (see Equation [3]). 



www.pu blish . csiro. au /journals/pasa 



11 



3xlO ie cm 



HO 



u=7 =Q r 25 



u=y=0A 



ti 

o 

-r-i 

m 

Oh 



LT =7=0,6 





-I — I 1 — I- 



I I I 



TT I 



I I I 



I I I 



3^ 



i i i 



n — i — i — i- 



-i— i- 




r 



■h- 1- 



V 



7.7xl0 16 cm 



-| 1 1 r 



I I I 



- HP I 

\ — I — h 



V 



I I I 



\ — I — h 




I I I 



I I I 



1— h-h 



V 




iK 



2xl0 17 cm 




Ik. 



0,1 



- - 0,01 



- -- 10 



-n-3 



10 



-4 



Figure 9: Correlations for Model B. The images along each row represent the 5 scenarios presented in 
Figure [5J namely a scattering region radius of 3xl0 16 cm with a separation distance of 6.03xl0 16 and 
1.66xl0 17 cm, 7.7xl0 16 cm with a separation distance of 2.00xl0 17 and 5.39xl0 17 cm, and 2xl0 17 cm with 
a separation distance of 3.70xl0 17 cm. Rows are then arranged in 3's, representing absorption at 77% 
(top), 91% (middle) and 97% (bottom). For each image the vertical axis shows polarisation while the 
horizontal axis shows magnification. 



12 



Publications of the Astronomical Society of Australia 







O 

P-i 



i so 
.9 i* 



0,01 10 3 10 4 

III I I — I [I II I I I I — I— 



I II I I I I l_ 



10 - 



3 4 6 a E 4 M £ 4 9 B 3 4GB M M 10 
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tiate between the two models. Note however that spec- 
tropolarimetric monitoring is required for high magni- 
fication events and an observational program to obtain 
high temporal sampling could be triggered based on 
simple photometric monitoring of the images. 

In addition, through the larger variations in polar- 
isation seen for all modelled scenarios of Model B, this 
study supports recent data for H1413+1143 (|Chae et al.l 
2001: [Chartas et aDl2004f ) which also suggests that a 
scattering region is most likely responsible for the in- 
creased polarisation found in BAL troughs. Such in- 
creased variation comes about in this model because 
the continuum region would be much more susceptible 
to microlensing than the larger scattering region. If 
this scenario is correct, then this hints that such a scat- 
tering region is also operating in other such quasars. 
For further work, the detailed temporal properties of 
the expected polarisation fluctuations need to be deter- 
mined with the goal of developing the most efficient ob- 
servational strategy to observe microlensed BAL quasars. 
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